Increasing the dimensionality of quantum walks using multiple walkers 
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We show that with the addition of multiple walkers, quantum walks on a line can be transformed 
into lattice graphs of higher dimension. Thus, multi-walker walks can simulate single-walker walks 
on higher dimensional graphs and vice versa. This exponential complexity opens up new applications 
for present-day quantum walk experiments. We discuss the applications of such higher-dimensional 
structures and how they relate to linear optics quantum computing. In particular we show that 
multi- walker quantum walks are equivalent to the BosonSampling model for linear optics quantum 
computation proposed by Aaronson & Arkhipov. With the addition of control over phase-defects 
in the lattice, which can be simulated with entangling gates, asymmetric lattice structures can be 
constructed which are universal for quantum computation. 



I. INTRODUCTION 

Quantum walks [TH3] (QWs), the quantum analogue 
of classical random walks, have emerged as an interest- 
ing approach to implementing quantum information pro- 
cessing tasks, and provide a useful tool for algorithm de- 
sign [5]. Recently there have been numerous experimen- 
tal demonstrations of optical QWs [6HT3]. Schreiber et 
al. pj [10] demonstrated a highly scalable optical QW on 
a line, whereby a single walker, simulated by weak coher- 
ent light, was temporally encoded and the size of the walk 
was limited only by loss. Shortly after, Broome et al. [8] 
demonstrated a single photon QW comprising an inter- 
ferometer constructed from wave-plates and birefringent 
crystals. A fully integrated device was demonstrated by 
Peruzzo et al. [9 in an interferometrically stable wave- 
guide. The state of the art is an experiment by Schreiber 
et al. [12 , who demonstrated a single walker quantum 
walk on a 2D lattice using temporally encoded position 
states. 

A central question in such demonstrations is how to 
scale the system to achieve the exponential complexity 
required for many interesting quantum information pro- 
cessing (QIP) [14] applications. This paper examines the 
ability to exponentially increase the complexity of such 
systems with the introduction of multiple walkers. The 
work by Peruzzo et al. [9] was a first step in experimen- 
tally achieving such complexity, and subsequently Refs. 
[TTJ [15] . However the pressing question is to understand 
how this can be employed for QIP tasks [13] [16j [17] . Un- 
derstanding such complexity is the central question we 
will address in this paper. 

Recently Rohde et al. [18 presented a formalism for 
multi- walker QWs. In this paper we extend this for- 
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malism and examine the applications for such higher- 
dimensional walks. In particular, we demonstrate that 
the introduction of multiple walkers is equivalent to em- 
ploying exponentially larger graph structures, one of the 
key ingredients in implementing QWs that cannot be ef- 
ficiently classically simulated, thereby opening the field 
to experiments of genuine interest to QIP. We demon- 
strate that multi-walker QWs are useful in simulating 
higher dimensional single- walker walks. Conversely, a 
single-walker walk on a higher dimensional graph struc- 
ture presents us with a versatile approach to mimicking 
multi- walker walks on simpler graphs. Such a simulation 
presents us with the tools necessary to understand multi- 
walker interactions, which are otherwise experimentally 
challenging to implement. This is particularly attractive 
to optical implementations, where directly implementing 
multi- walker interactions is difficult, requiring massive 
optical non-linearities or entangling gates [19j. In light 
of the experiments by Peruzzo et al. [9 and Matthews 
et al. [11 we also give consideration to entangled input 
states and discuss how they may, or may not, benefit 
experiments. 

There are two general classes of quantum walks that 
have been shown to be universal for quantum computa- 
tion: continuous- (by Childs [5]) and discrete-time (by 
Lovett et al. [20]) QWs. In this paper we will focus 
on discrete- time QWs. We demonstrate that multiple 
walkers generate a virtual graph, whose dimensionality 
is higher than that of the underlying graph. For exam- 
ple, two walkers transforms a linear graph into a regular 
lattice in two dimensions, and with three walkers into a 
three-dimensional lattice. Thus, as the number of walkers 
increases, the size of the virtual graph increases exponen- 
tially. This is an interesting observation, as it implies that 
with frugal physical resources we can construct a system 
of exponential complexity. 

In particular, we demonstrate that a multi- walker QW 
is isomorphic to the recently introduced BosonSam- 
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pling model for optical quantum computation (QC) in- 
troduced by Aaronson & Arkhipov [21]. This model is 
strongly believed to not be universal for QC, but is 
nonetheless believed to be classically hard to simulate, 
making it of direct interest to optical QC experimental- 
ists. However, when lattice defects are introduced, which 
can be simulated with entangling gates, the scheme be- 
comes universal. Furthermore, we argue that a multi- 
walker QW on an efficiently sized graph (equivalently 
BosonSampling) is universal for QC only if a QW on 
an exponentially large graph with N vertices is universal 
with control over just 0(log N) coins, which would be a 
surprising result, giving further weight to the argument 
that BosonSampling is not universal for QC. However, 
efficient QC is possible with the addition of CPHASE 
gates, which will be discussed. We also discuss the poten- 
tial to using entangled state preparation to overcome the 
difficulty of implementing entangling operations within a 
QW. 



II. THE QUANTUM WALK FORMALISM 

A QW consists of a walker, which may occupy some 
number of position states in a graph. In a discrete-time 
QW, a walker is a bipartite system consisting of a posi- 
tion (x) and a coin (c) value, |x,c). The position literally 
corresponds to the walker's position in the graph, i.e. at 
which vertex the walker resides, whereas the coin is an 
ancillary parameter whose value determines the direction 
of propagation of the walker. The walk proceeds using 
two operations - coin (C) and step (S). C coherently 
randomises the value of the coin parameter, whereas S 
uses the value of the coin to propagate the walker. In a 
classical random walk the evolution proceeds by flipping 
a classical coin and then depending on its value propa- 
gating the walker in the appropriate direction. In a QW, 
on the other hand, the coin may be in a coherent super- 
position. Thus after propagation the walker is, in general, 
in a superposition of different positions. 

As an example, we consider the archetypal example of 
a quantum walk - a single- walker on a linear graph. Then, 
using a so-called Hadamard coin, the evolution proceeds 



III. THE MULTI- WALKER QUANTUM WALK 
FORMALISM 



Following Ref. [18 , we represent multiple walkers on 
arbitrary graph structures using walker operators. For- 
mally, we define a walker operator of the form w(x,c)^ . 
Previous authors have given consideration to the differ- 
ences between fermionic and bosonic Q Ws [TTJ [13] . While 
such differences are valuable to understand, in this paper 
we will focus on bosonic QWs, given their applicability 
to present-day optical implementations. Thus, w(x,c)^ 
can be directly interpreted as a bosonic creation opera- 
tor. The same formalism could be applied to fermionic 
walkers, but with different commutation relations. 

We define the coin and step operators as, 

C(t): w(x,c)^ ^A ( *\t)w(x,j)\ 

S : w(x, j) ] ^ w(j,x) ] . (2) 

Here n x is the neighbourhood of position x, i.e. the ver- 
tices to which x is connected in the graph, and is 
an operator which determines the evolution of the coin 
value within the respective neighbourhood. An example 
is illustrated in Fig. Ill Note that w(x, c)t|0) forms a basis 
for all single- walker states. Here |0) represents an empty 
graph, where no walkers are present. In an optical con- 
text this corresponds to the vacuum state. The evolution 
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FIG. 1: Example graph with 5 positions. A unitary coin op- 
erator acts on the neighbourhood of each vertex i. 

See also Ref. [22 for a discussion on multi-particle 
quantum walks. 



IV. VIRTUAL WALKERS 



c\x,±i) = (M)±k,-i>)/v% 

S\x,c) = \x + c,c). 



(i) 



Evidently, at each step, the walker enters a superposition 
of the position states immediately to its left and right. 
The QW has some interesting features compared to clas- 
sical walks, including localisation and an enhanced rate- 
of-spread of the walker [3 j. The increased rate of spread 
of a QW compared to the classical counterpart is one of 
the features of QWs that makes it of direct interest to 
certain algorithmic applications such a search algorithms. 



With n walkers in the system the basis states will be of 
the form w(x\, ci)^w(x2, c 2 y . . . w(x n , c n )^|0). We conse- 
quently define virtual walker operators of the form, 

v(x, c) ] = w(x 1 ,c 1 ) ] w(x 2 , c 2 ) f . . . w(x n , c n ) ] (3) 

(up to normalisation) where x — {#i, x 2l • • • •> x n} an d c = 
{ci, C2, . . . , c n }. We refer to these as virtual walkers since 
each set of walker creation operators, while represent- 
ing multiple particles, can be thought of as a single ex- 
citation, inheriting their properties from the underlying 
bosonic creation operators - they are indistinguishable, 
commutative and obey (0\v(x,c)v(y,dy\0) = 
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Let us consider a linear graph as shown in Fig. [2^ top) . 
With two walkers in the system at positions (i, j), we can 
immediately read off which transitions between virtual 
walkers are allowed, Fig. [2|bottom). With two walkers it 
can be seen that the pair of walkers at vertices (i, j) may 
transition to (i + 1, j + 1), + — 1), (i — 1, j + 1) 
and (i — 1, j — 1). For a linear graph with identical coins 
at all positions, all are equal. However, in general 
they may be distinct, allowing arbitrary graphs struc- 
tures to be defined with coins of different dimension. 

From the allowed transitions we can construct a vir- 
tual graph. In the example of an underlying linear graph, 
Fig. |3|top), it can be seen that the virtual graph is a 
two-dimensional lattice graph, Fig.[3jbottom). Thus, the 
virtual graph, with a single walker, exhibits the same al- 
lowed transitions as the underlying graph with two walk- 
ers. Generally, in the virtual graph there will be an edge 
(i, j) — » (&,/) if in the underlying graph there are edges 
i — » k and j —> I. Representing the virtual walkers in a 
two-dimensional space we can see that the graph is com- 
posed of overlapping four pointed stars as the basic build- 
ing block. Thus, with the addition of a second walker, a 
linear walk consisting of one walker is transformed into a 
two-dimensional lattice graph consisting of a single vir- 
tual walker. With the addition of a third walker the ba- 
sic building block will be a three-dimensional star within 
a three-dimensional lattice, and for n walkers the vir- 
tual graph will be an n-dimensional lattice. As the num- 
ber of walkers increases, clearly the size of the virtual 
graph increases exponentially, even though the underly- 
ing linear graph is not increasing in size. Thus, the vir- 
tual graph allows us to trade quantum complexity for 
exponentially growing classical complexity. We have fo- 
cussed on linear graphs as a simple example, since most 
present-day optical implementations have been restricted 
to linear graphs. The described walker operator formal- 
ism may be applied to studying arbitrary graph struc- 
tures. More generically, a connected graph K\ G \ on graph 
G with n walkers, will transform into the connected graph 
K\ G \ n . The de gree of vertex p Xl , X2 ,..., Xn : ° r equivalently 
the size of the coin space, in the virtual graph is related 
to the degree of the vertices q y in the underlying multi- 
walker graph by \p XliX2i ... lXn \ = EIILi wmcn in gen- 
eral grows exponentially with n. 
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i - 1 i + lj-1 j + l 

FIG. 2: Evolution of two walkers at time t to time t + 1. Black 
positions are potentially occupied, while white positions are 
unoccupied. We assume the two walkers are in positions (i, j). 
Then the allowed transitions can be immediately read off the 
graph. 

In Fig. [4] we show the walker probability plots for two 




FIG. 3: (top) A linear walk with 6 vertices. Numbers denote 
different positions, (bottom) The equivalent walk where two 
walkers are mapped to a single virtual walker. Numbers rep- 
resent pairs of positions. Note that the graph is symmetric 
about the diagonal axis. This is because the particles are 
indistinguishable, therefore (xi,X2) = (x2,xi), and all non- 
diagonal vertices are 2-fold redundant. 

equivalent situations: (1) P(x\,X2) is the position of a 
single walker on a 2D lattice of the form shown in Fig. 
|3j where x\ and X2 represent the two spatial dimensions; 
and, (2) P(x\,X2) is the coincidence probability of two 
walkers on a linear graph, where x\ and x^ denote two 
distinct positions on a line. Note that because of the 
isomorphism, virtual walkers on a line exhibit the same 
properties as a single walker on a lattice. In the graphs 
in Fig. [4j ballistic spreading is observed, an archetypal 
property of quantum walks [3]. 

The idea of mapping an efficiently sized graph to an ex- 
ponentially large one was very recently described by Un- 
derwood & Feder [23] , in which they showed that such 
such systems can be made universal for quantum com- 
putation using Bose-Hubbard type interactions between 
interacting bosons. We will later demonstrate a similar 
result for the universality of our virtual graphs via the 
introduction of defects into the virtual graph. 

V. INTRODUCING LATTICE DEFECTS 

We now consider some interesting properties of this 
mapping. First, note that in an arbitrary linear optics 
network a single photon input state will give the same 
single click-statistics as weak coherent light. Thus, while 
two walkers on a line must be implemented using two 
photons, one walker on a lattice can be implemented with 
either a single photon or with coherent light. Note that 
the latter system is purely classical and the former is, 
in general, highly entangled. Thus, a 2D spatial lattice 
represents a system, which, while containing no entan- 



FIG. 4: (Color online) Probability plots for one walker in a lattice of two spatial dimensions, or equivalent ly, the coincidences 
for two walkers on a line (a virtual walker). Because the line is translation invariant, all coins are equal and are defined 
as Hadamard operations. Note the walkers /virtual walker exhibit ballistic spreading, a unique feature of quantum mechanical 
walks. 



glement, is able to simulate a system with entanglement. 
Importantly however, even though no entanglement is 
present in the latter system, coherence is still required. 

Next one might ask the question whether entangling 
operations can be implemented in such a system. Suppose 
we wish to implement a generalisation of the CPHASE 
gate, which applies a phase- flip at position {xi, . . . , x n } 
in the virtual graph only if there is a walker at all 
positions X\,...,x n in the underlying n- walker graph. 
For example, a two-particle CPHASE gate would be 
defined as CPHASE^ jClj;E2 ,c2 : w(xi, Ci)^w(x2, £2)^ — ^ 
—w(xi,CiYw(x2 : C2)t. In the linear graph with two walk- 
ers, this clearly requires an entangling operation - in fact, 
in an optical context, it would require a massive non- 
linear interaction, putting a ir phase-shift on two-walker 
terms. However, in the 2D lattice with a single walker this 
is no longer the case. Rather, we can pick out an individ- 
ual term of interest and apply a non-entangling PHASE 
gate at the respective position in the virtual graph (see 
Fig. [5|. We refer to applications of such localised PHASE 
gates as lattice defects. Such defects present us with sev- 
eral interesting properties: 

1) We can simulate entanglement with classical states, 
by accepting an increase in the graph dimensionality. 
This is interesting from a quantum simulation perspec- 
tive. Specifically, by accepting an exponential overhead 
(against the number of walkers) in graph size, we can 
simulate a multi-walker QW with entangling operations, 
such as CPHASE gates and non-linear interactions, even 
though no entanglement is present. 

2) We can use CPHASE gates on a graph with mul- 
tiple walkers to efficiently 'etch' an arbitrary pattern of 
phase-defects into the corresponding virtual graph. This 
is particularly interesting as it allows us to break the 
translation symmetry otherwise present in the virtual 
graph, opening the possibility to explore more elabo- 
rate, non-trivial graph structures with potentially more 
interesting dynamics. The virtual graph will always be 
symmetric about the diagonal axis, since with indistin- 
guishable walkers {^i,^} = {^2>^i}- However, modulo 
this constraint, arbitrary structures of phase-defects can 



be efficiently implemented with an appropriate choice of 
CPHASE gates. 

3) Vertices on the main diagonal of the virtual graph 
play a special role, since these correspond to both walkers 
being in the same position. Thus, a phase shift applied 
to these vertices represents a non-linear interaction where 
the applied phase is a function of how many walkers re- 
side at the vertex. 



jCPHASEj \ I j, 




FIG. 5: An example mapping between entangling CPHASE 
gates on a linear graph with two walkers and non-entangling 
PHASE gates (denoted by '-') in the corresponding virtual 
graph. An appropriate choice of CPHASE gates allows us to 
efficiently 'etch' arbitrary phase-defects into the correspond- 
ing virtual graph, subject to the constraint that the virtual 
graph is symmetric about the diagonal axis. This allows us 
to generate a larger class of virtual graphs that are in general 
not translation-symmetric and which may exhibit less trivial 
behaviour. 

4) With the addition of generalised CPHASE gates 
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FIG. 6: (Color online), (left) Mapping an N = 3 quantum walk to a beamsplitter network, (right) Mapping an arbitrary single 
beamsplitter operation between modes (2, 2) and (3, 2), of 9 modes, to an N = 3 quantum walk. The routing of the modes to 
which the beamsplitter operation is applied are shown in bold red. First, a choice of permutation coin operators maps the two 
modes to have the same coin value, after which the step operator maps them to the same position. Then an additional coin 
applies an arbitrary beamsplitter operation to the two modes. A subsequent step and coin route the two modes back to their 
original position/coin values. 



the model becomes universal for QC. To see this we 
demonstrate that the model can be directly mapped to a 
qubit system with a universal gate set. Let us begin with 
a linear walk. Then at each vertex there are two allowed 
coin states (left and right), w(x, x + 1)1" and w(x, x — 1)^. 
Thus, with a single walker at each vertex in a superpo- 
sition of the two allowed coin states, we define a set of 
qubits - one at each vertex, a-w(x, x-\-l)^ +j3-w(x, x— 
Arbitrary single qubit operations can trivially be per- 
formed with an appropriate map within the subsequent 
coin operations (C). The final step is to introduce the 
CPHASE gate, the last remaining required gate to form 
a universal gate set. To implement this gate we begin 
with a 'routing' argument. As an example, Fig. fright) 
illustrates how to map two arbitrary wires to lie within 
a single coin operator C. Now, we allow this C to be 
more than just a unitary map on the creation operators, 
but a CPHASE operator acting on the appropriate rails 
within the 'bundle' of modes corresponding to a given po- 
sition (Fig. |6fleft)). Note that implementing CPHASE 
gates directly is in general non-deterministic when imple- 
mented using just linear optics [19]. Thus, this final step 
remains the challenging one. 

Interestingly, this implies that, when looked at from 
the perspective of the corresponding virtual graph, an 
exponentially large regular lattice graph is universal for 
QC, with an appropriate choice of phase-defects. Without 
phase-defects the model is unlikely universal for QC since 
it does not allow for the construction of arbitrary graphs 
and therefore cannot be mapped to a protocol such as 
that by Lovett et al [2D] . This result on the other hand 
suggests that full control over the graph is not necessary. 
Rather, a regular lattice is sufficient, with full control 
over phase-defects. 

In the same way that the QW formalism provided a 
new approach to visualising and designing quantum al- 
gorithms in a graph theoretic context, this alternate ap- 
proach, based on regular lattices with phase-defects, may 
provide insight into new approaches to algorithm design. 



In Ref. [12] Schreiber et al. make use of the aforemen- 
tioned isomorphism to simulate two walkers on a line. 
They begin by experimentally demonstrating an optical 
2D QW, and with the addition of an electro-optic mod- 
ulator (EOM) are able to apply arbitrary phase defects 
into the 2D lattice. They subsequently use the isomor- 
phism to simulate two walkers, and are able to simulate 
both entangling and non-linear interactions, providing a 
rich test-bed for the simulation of multi-particle dynam- 
ics. With appropriate choices of coins and defects they are 
able to simulate highly entangled two-particle dynamics 
using just classical light. 



VI. IMPLICATIONS FROM COMPLEXITY 
THEORY 

We now relate the described multi- walker formalism to 
two protocols and comment on the implications for com- 
plexity theory. Having discussed the potential for multi- 
walker QWs to achieve exponential complexity, the ques- 
tion is whether this isomorphism can be utilised to per- 
form universal QC, or at least some interesting problems 
which are intractable on a classical computer. We relate 
the QW formalism to the protocols by Knill, Laflamme 
& Milburn [19], and Aaronson & Arkhipov [21 . With the 
addition of feedforward the multi-walker QW formalism 
can be made universal for QC, and in the absence of 
feedforward is equivalent to BosonSampling. 



A. Universal quantum computation using quantum 
walks? 

Childs [5 first demonstrated that a continuous time 
QW is universal for QC, and a similar proof exists in 
the discrete time case [20]. Here an exponentially sized 
graph with 0(2 n ) vertices is required to simulate a qubit 
circuit comprising n qubits. Of course, in a physical sys- 
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tern where position states literally correspond to a phys- 
ical position, an exponential number of vertices is clearly 
not possible [24 j. However, we have demonstrated that 
(even without phase-defects) with the addition of mul- 
tiple walkers the complexity of the system grows expo- 
nentially, in terms of the number of available basis states 
(see Ref. [21] for further discussion on the complexity of 
linear optics systems with multiple photons), raising the 
question as to whether this exponential complexity can 
be mapped to the Childs [5] or Lovett et al. [20 approach 
and is universal for QC. 

In the work by Childs [5] and Lovett et al. [20], ar- 
bitrary graph structures are required for universal QC. 
There, widgets implementing elementary gate operations 
are attached to 0(2 n ) quantum wires, which propagate 
basis states. For example, there is a widget which imple- 
ments a CNOT gate, and widgets for implementing the 
remaining single qubit operations necessary for universal 
QC. 

Because the number of vertices in the virtual graph is 
exponentially larger than the underlying graph, and we 
only have direct control over the underlying graph, it is 
clear that not all combinations of virtual graphs may be 
constructed. We do not have sufficient control over the 
isomorphism to generate arbitrary virtual graphs. Ad- 
ditionally, to fully utilise an exponentially large virtual 
graph one must have control over an exponential num- 
ber of coin parameters. Since we can only directly im- 
plement coins on the underlying graph it is clear that 
not all combinations of coins in the virtual graph can be 
implemented. We therefore suggest that, without phase- 
defects, universal QC is unlikely using the mult i- walker 
formalism since it cannot be mapped to arbitrary ex- 
ponentially large single-walker walks. Indeed, a simple 
counting argument gives stronger weight to the claim 
that multi-walker QWs cannot replicate a Lovett et al. 
[20] protocol. Suppose the underlying graph contains k 
vertices. Then, there are k locations at which coins can 
be applied. Similarly, in the corresponding virtual graph 
there are also k allowed unique combinations of locations 
where coins may be applied. In the Lovett et al. [20] ap- 
proach, the walk is defined by an exponential number 
of locations at which coins may be applied. Thus, if all 
exponentially sized graphs are to be simulated then, if 
they can be simulated at all, an exponential number of 
time steps will be required simply to allow for the desired 
combinations of coin operators, suggesting that efficient 
universal QC is not possible. If efficient universal QC is 
possible on a virtual lattice without phase defects, this 
would imply that a quantum walk with N basis states is 
universal with control over just 0(log N) coins. 

Note that the above discussion is not a formal proof 
that a multi-walker quantum walk is not universal for 
quantum computing, but rather an intuitive plausibility 
argument based on a counting approach. 



B. Relation to linear optics quantum computing 

Although unlikely to be universal for QC (e.g. as 
proposed by Knill, Laflamme & Milburn (KLM) [13]). 
mult i- walker QWs can efficiently implement the Boson- 
Sampling algorithm proposed by Aaronson & Arkhipov 
(AA) [2T] , We first show that multi- walker QWs can ef- 
ficiently simulate arbitrary linear optics networks, and 
then relate this to KLM and AA. 



1. Mapping quantum walks to optical networks 

From Eq. [2] it is easy to write down a mapping of a 
quantum walk to a beamsplitter network. Note that the 
coin operator acts individually on each 'bundle', where 
a bundle of modes represents a walker position. Reck et 
al. [25] showed that an arbitrary unitary map of the form 
a \ — > . UijOj has an efficient decomposition into beam- 
splitters and phase-shifters. Thus the operation act- 
ing on each bundle can be efficiently constructed. Then it 
can easily be seen that Fig.^left) is an efficient decompo- 
sition of an arbitrary quantum walk into a beamsplitter 
and phase-shifter network. 

2. Mapping arbitrary optical networks to a quantum walk 

Mapping in the opposite direction, from a beamsplitter 
and phase-shifter network to a quantum walk is not as 
trivial. We aim to show that a QW can act as a decom- 
position for an arbitrary linear optics network. It suffices 
to show that an efficient number of coin and step op- 
erators can be used to apply beamsplitter/phase-shifter 
operations between arbitrary pairs of modes. To do this 
we first use permutations within the first coin opera- 
tor. These permutations, followed by the step operator, 
route two arbitrary modes to the same position 'bundle'. 
Once both modes have been mapped to the same po- 
sition an arbitrary beamsplitter /phase-shifter operation 
can be applied to the corresponding two modes with an 
appropriate choice of Then, another application of 
the same step and permutation operators route the re- 
spective modes back to their original coin and position 
states. In general, multiple operations can be performed 
in parallel. Formally, the decomposition can be expressed, 

-B X iCi,X2C2 = Px\l,X\C\ ' -Px2l,X2C2 ' ^ ' 

B\ X \,lX2 ' ^ ' Px\l,X\C\ ' Px2l,X2C21 (^) 

where Pij y u is a permutation between modes (i, j) and 
(fc, I), S is the usual step operator, and Bij^i is the beam- 
splitter operation between modes and (fc, I). Here P 
and B are the choices of coin operators, so the sequence 
of operators still consists of iterations of coin and step 
operators. An example is illustrated in Fig. fright). Fi- 
nally, we know from Reck et al. [25] that a polynomial 
number of iterations of coin and step operators are now 



7 



sufficient to decompose an arbitrary linear optics network 
into a QW. Thus, a QW is an efficient decomposition for 
arbitrary unitary maps. 



3. Relation to KLM and AA 

We now briefly discuss the relationship between multi- 
walker QWs and LOQC. We discuss two different ap- 
proaches to LOQC - KLM and AA. These two models 
are subtly different. The KLM approach requires arbi- 
trary linear optics networks and the addition of multi- 
photon input states, measurement, post-selection and 
fast feed- forward. This approach is known to be univer- 
sal for QC (i.e. it can solve BQP-complete problems). 
The AA approach on the other hand also requires gen- 
eral linear optics networks and multi-photon input states, 
but mitigates the necessity for post-selection and feed- 
forward. This approach implements the so-called Boson- 
Sampling problem, which is strongly believed to be 
classically hard (residing in the class BosonSampP- 
complete) [21] . While far simpler and significantly exper- 
imentally less demanding, this approach is neither known 
nor believed to be universal for QC [21]. Note that both 
approaches require arbitrary beamsplitter networks. We 
have shown in the previous section that QWs act as a de- 
composition for arbitrary linear optics networks. Thus, it 
follows that multi- walker QWs can be used to implement 
both KLM and AA. 

In the case of KLM the QW must be complemented by 
measurement and feed-forward. Thus the corresponding 
QW consists not only of coin and step operators, but also 
measurement operators, and the evolution of the system 
can be represented as ]\ t M(t) • 5 • C(£, m t -i), where M 
is a measurement operator, and m t denotes the measure- 
ment outcome at time t. 

AA on the other hand can be implemented using a QW 
in the more conventional way of iterative applications of 
just coin and step operators. We have demonstrated that 
in the absence of feedforward, multi-walker QWs are un- 
likely to be universal for QC, since there is no general 
mapping to the QW protocol of Lovett et al [20 . This 
gives further weight to the claim that BosonSampling 
is likely not universal for QC. As discussed earlier, if 
BosonSampling were to be universal for QC, this would 
imply that the protocol of Lovett et al. [20] is universal for 
QC even when acting on exponentially sized regular lat- 
tice structures of size N with control over only O (log TV) 
coins. We can establish the condition that BosonSam- 
pling is universal for efficient QC only if a single- walker 
walk on exponentially large graphs can be implemented 
using time-dependent periodic graph structures, with an 
efficient number of time steps. Conversely, because of the 
isomorphism between multi-walker QWs and arbitrary 
linear optics networks, the work of A A provides a useful 
application for multi- walkers QWs. Although such regu- 
lar lattice structures are unlikely universal for QC, the 
question is raised whether nonetheless some interesting 



problems can be solved on such structures. 



VII. ENTANGLED INPUT STATES 

In conventional approaches to QC, such as the circuit 
and cluster state [26 , 27 models, it is known that initially 
entangled states act as a useful resource to the benefit of 
QC. Also they provide an advantage in QW algorithms, 
as for the example the graph isomorphism problem [28J. 
One might ask the question whether the same applies 
in the context of our multi- walker QW model. Specifi- 
cally, we know that with the addition of C PHASE gates 
our model becomes universal. However, such entangling 
gates are both challenging and non-deterministic in the 
optical context. Thus, the question is whether the non- 
trivial gates can be commuted to the beginning of the 
circuit such that the challenging part of the problem is 
transferred to a resource state preparation problem after 
which the algorithm is 'easily' implemented using just 
the passive, and deterministic, sequence of coin and step 
operators. Here we show that this is not the case. We 
present a simple example of a circuit which illustrates 
that, in general, such resource state preparation is as dif- 
ficult as implementing the CPHASE gates within the 
circuit. Thus, the idea of commuting CPHASE gates to 
the beginning of the circuit and treating it as a resource 
state preparation problem is of little practical benefit. 

In Fig. [7] we present an extremely simple circuit con- 
sisting of two control wires and two additional wires to 
which two consecutive beamsplitters are applied. After 
each beamsplitter comes a CPHASE gate, controlled by 
one of the two control wires. Using the commutation re- 
lations in Fig. [7jinset top) we commute the CPHASE 
gates to the beginning of the circuit and treat all the com- 
ponents preceding the QW as a resource state | ^resource)- 
Evidently, the number of controlled gates in the state 
preparation procedure is as large as the number of con- 
trolled gates we had in the original circuit. On one hand 
this suggests that resource state preparation does not 
mitigate the difficulties inherent in implementing con- 
trolled gates within the circuit. On the other hand, be- 
cause all the non-deterministic gates occur at the begin- 
ning of the circuit, they can be implemented 'offline', a 
trick which is well known for optical cluster states [29] . 
However, in the optical context, entangling controlled 
gates are extremely challenging. Thus we suggest that 
the trick of resource state preparation does not represent 
an efficient solution to the problem. 

It should also be noted that while in our simple toy ex- 
ample commuting all controlled gates to the beginning of 
the circuit is possible, in general this is not the case. It is 
trivial to draw circuits where the described commutation 
relations do not allow for all entangling operations to be 
commuted to the beginning of the circuit. For example, 
consider the simple circuit shown in Fig. [8] Here no trivial 
commutation relation exists enabling the CPHASE gate 
to be commuted to before the beamsplitters. Commu- 
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tation would involve a four-mode entangling operation. 
Nonetheless, it may be the case that there is a subset of 
interesting problems for which commutation is possible 
and the resource state preparation is more practical than 
implementing controlled gates within the circuit. 
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FIG. 7: (inset top) Commutation relation for a 7r-phase gate 
and a beamsplitter, and a swap gate and a beamsplitter, (bot- 
tom) Commuting controlled-7r gates within a trivial QW to 
the beginning of the circuit, such that the circuit consists of 
an entangled resource state and a QW without any controlled 
gates. 



tially growing virtual lattice. Although the complexity of 
the system grows exponentially, it is unlikely to be uni- 
versal for quantum computation as the resulting graph 
has inherent symmetries that limit our ability to engineer 
arbitrary graph structures, unless CPHASE gates are in- 
troduced, which generate arbitrary lattice defects. How- 
ever, in the absence of CPHASE gates, a multi-walker 
quantum walk is isomorphic to the BosonSampling 
problem, which is of practical interest, and within reach 
of near-term experiments. Our arguments give further 
weight to the assertion that BosonSampling is not uni- 
versal for quantum computation, since there is no equiv- 
alence to the scheme for universal quantum computation 
using quantum walks. Nonetheless, multi- walker QWs are 
likely capable of tasks classically intractable and there- 
fore of practical interest. We demonstrated that single- 
walker 2D walks can simulate two- walker ID walks and 
therefore a strictly non-entangling system (i.e. one walker 
simulated with weak coherent light) can simulate a highly 
entangling system. This observation paves the way for 
elementary demonstrations of quantum walks with mul- 
tiple walkers, via the addition of extra dimensions. We 
also demonstrated that with the addition of entangling 
gates, such as CPHASE gates, the symmetry of the cor- 
responding virtual graph can be broken by adding defects 
to the lattice, paving the way for non-trivial graph struc- 
tures to be studied. Of particular interest, with the addi- 
tion of CPHASE gates, or equivalent ly phase-defects in 
the corresponding virtual graph, the system becomes uni- 
versal for quantum computation, despite the fact that we 
do not have full control over the graph structure. Rather, 
a regular lattice combined with full control over phase- 
defects is sufficient for universal quantum computation. 



X 



/ — 1 


1 




b— 



FIG. 8: A simple example of a beamsplitter network with a 
single CPHASE gate between two modes. There is no triv- 
ial commutation relation to commute the CPHASE gate to 
before the beamsplitters. Commutation would require a four- 
mode entangling operation. 



VIII. CONCLUSION 

We have shown that adding walkers to a quantum walk 
system is equivalent to a single walker on an exponen- 



Acknowledgments 



We thank Timothy Ralph, Scott Aaronson and Aurel 
Gabris for helpful discussions. This research was con- 
ducted by the Australian Research Council Centre of Ex- 
cellence for Engineered Quantum Systems (Project num- 
ber CE1 10001013). We acknowledge support by the grant 
MSM 6840770039 of the Czech Ministry of Education and 
the Doppler Institute, FNSPE CTU in Prague. 



[1] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. [2] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, 
A 48, 1687 (1993). 



9 



in Proceedings of the 33th STOC, New York p. 50 (2001). 
[3] J. Kempe, Cont. Phys. 44, 307 (2003). 
[4] S. E. Venegas-Andraca (2012), arXiv:1201.4780. 
[5] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). 
[6] H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Moran- 

dotti, and Y. Silberberg, Phys. Rev. Lett. 100, 170506 

(2008). 

[7] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 

P. J. Mosley, E. Andersson, I. Jex, and C. Silberhorn, 

Phys. Rev. Lett. 104, 050502 (2010). 
[8] M. A. Broome, A. Fedrizzi, B. P. Lanyon, I. Kassal, 

A. Aspuru-Guzik, and A. G. White, Phys. Rev. Lett. 

104, 153602 (2010). 
[9] A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, 

A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, 

K. Worhoff, et al., Science 329, 1500 (2010). 
[10] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 

I. Jex, and C. Silberhorn, Phys. Rev. Lett. 106, 180403 

(2011). 

[11] J. C. F. Matthews, K. Poulios, J. D. A. Meinecke, 
A. Politi, A. Peruzzo, N. Ismail, K. WorhofT, M. G. 
Thompson, and J. L. O'Brien (2011), arXiv:1106.1166. 

[12] A. Schreiber, A. Gabris, P. P. Rohde, K. Laiho, 
M. Stefanak, V. Potocek, I. Jex, and C. Silberhorn, Sci- 
ence 336, 55 (2012). 

[13] L. Sansoni, F. Sciarrino, G. Vallone, P. Mataloni, 
A. Crespi, R. Ramponi, and R. Osellame, Phys. Rev. 
Lett. 108, 010502 (2012). 

[14] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 



Cambridge, 2000). 
[15] J. O. Owens, M. A. Broome, D. N. Biggerstaff, M. E. 
Goggin, A. Fedrizzi, T. Linjordet, M. Ams, G. D. Mar- 
shall, J. Twamley, M. J. Withford, et al., New J. Phys. 
13 (2011). 

[16] A. Aspuru-Guzik and O. Walther, Nature 8, 285 (2012). 
[17] S. D. Berry, New J. Phys. 13, 033037 (2011). 
[18] P. P. Rohde, A. Schreiber, M. Stefanak, I. Jex, and C. Sil- 
berhorn, New J. Phys. 13, 013001 (2011). 
[19] E. Knill, R. Laflamme, and G. Milburn, Nature (London) 

409, 46 (2001). 
[20] N. B. Lovett, S. Cooper, M. Everitt, M. Trevers, and 

V. Kendon, Phys. Rev. A 81, 042330 (2010). 
[21] S. Aaronson and A. Arkhipov, STOC (2011), 

arXiv:1011.3245. 
[22] K. Mayer, M. C. Tichy, F. Mintert, T. Konrad, and 

A. Buchleitner, Phys. Rev. A 83, 062307 (2011). 
[23] M. S. Underwood and D. L. Feder (2012), 

arXiv: 1204.6021. 
[24] P. P. Rohde (2010), arXiv:1010.4608. 
[25] M. Reck, A. Zeilinger, H. J. Bernstein, and P. Bertani, 

Phys. Rev. Lett. 73, 58 (1994). 
[26] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 

5188 (2001). 

[27] R. Raussendorf, D. E. Browne, and H. J. Briegel, Phys. 

Rev. A 68, 022312 (2003). 
[28] S. D. Berry and J. B. Wang, Phys. Rev. A 83, 042317 

(2011). 

[29] M. A. Nielsen, Phys. Rev. Lett. 93, 040503 (2004). 



